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ABSTRACT 

Hamilton et al. have suggested an invaluable scaling formula which describes how the 
power spectra of density fluctuations evolve into the nonlinear regime of hierarchical 
clustering. This paper presents an extension of their method to low-density universes 
and universes with nonzero cosmological constant. We pay particular attention to 
models with large negative spectral indices, and give a spectrum-dependent fitting 
formula which is of significantly improved accuracy by comparison with an earlier 
version of this work. The tendency of nonlinear effects to increase power on small scales 
is stronger for spectra with more negative spectral indices, and for lower densities. 
However, for low-density models with a cosmological constant, the nonlinear effects 
are less strong than for an open universe of the same f2. 
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1 INTRODUCTION 

The power spectrum of density fluctuations is a statistic of 
central importance in cosmology, as it describes a combi- 
nation of the primordial deviations from homogeneity and 
their subsequent modification by the matter content of the 
universe. One practical obstacle to a measurement of this 
useful function is that the present universe occupies a state 
where nonlinear gravitational growth of density fluctuations 
has altered the form of the fluctuation spectrum. 

Prior to 1991, it was assumed that this problem was 
tractable only in two extreme limits. On large scales, linear 
theory applies and the power spectrum scales as the square 
of the density growth factor. For very large wavenumbers, we 
enter the 'stable clustering' regime, where there is a simple 
relation between the power-law index of the linear spectrum 
and that of the output nonlinear spectrum (see Section 73 
of Peebles 1980) . It was therefore a significant breakthrough 
when Hamilton et al. (1991; HKLM) suggested a scaling pro- 
cedure which allowed an accurate description of the transi- 
tion regime between these two limits in terms of an empirical 
universal function. 

In a previous paper (Peacock & Dodds 1994; PD), we 
extended the HKLM procedure in a number of ways. First, 
we presented a version of the method which worked with 
power spectra, rather than HKLM's choice of integrated cor- 
relation function. Second, we considered the modifications to 
the HKLM argument needed to work in a universe of arbi- 
trary density, rather than HKLM's fi = 1. We applied this 
extended HKLM method to a compilation of clustering data, 
and concluded that the linearized data for k < 0.4/iMpc -1 
(h = -ffo/100 kms -1 Mpc -1 ) were consistent with a spectral 
index of about —1.5 on small scales, steepening to close to 
the primordial n = 1 on large scales (-P(fe) cx k n ). 

However, the universal scaling procedure of HKLM was 
known to fail for more negative spectral indices, as has been 



investigated in some detail recently for Q = 1 models by 
Jain, Mo & White (1995; JMW). A common description of 
the linear power spectrum is in terms of the CDM model, 
whose power-law slope tends to —3 on very small scales, so 
the HKLM method and its extensions in PD may not work 
very well if we attempt to investigate scales significantly 
smaller than those studied in PD. We have therefore run 
an additional ensemble of A-body simulations which con- 
centrate on the case of spectra with n < — 1 and the CDM 
spectrum. It turns out that a simple slope-dependent cor- 
rection can be made to the PD formulae which provides an 
excellent description of the nonlinear data over essentially 
all regimes of interest, and this is described below. 



2 THE HKLM METHOD 

The key argument of HKLM is that gravitational collapse 
causes a change of scale. By regarding the integrated cor- 
relation function £(r) as measuring the number of excess 
neighbours within radius r, they suggested that observed 
nonlinear correlations on nonlinear scale r NL be related to a 
pre-collapse linear scale via 



(1) 



HKLM then conjectured that, having translated scale, the 
linear and nonlinear correlations had some universal relation 



Cnl^nl) = /nl[£l(7l)]. 



(2) 



The function / NL must behave as / NL (x) = x in the linear 
i«l limit and / NL (x) oc x 3 ^ 2 in the stable-clustering x ^> 1 
limit, and must be determined numerically around x ~ 1. 

PD argued that a very similar argument could be made 
to work for power spectra, using a dimensionless version of 
the power spectrum: A 2 is the contribution to the fractional 
density variance per unit In k. In the convention of Peebles 
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(1980), this is 
A 2 (fc) = da " 



V 



dlnfc (2tt) : 



■47rfc 3 |<5 fe 



(3) 



(V being a normalization volume), and the relation to the 
correlation function is 



2 dfc sinfcr 
fe kr 



(4) 



Since £(r) can be thought of as A 2 at some effective 
wavenumber, this suggests the fc-space version of HKLM: 

fc L = [1 + A 2 , L (fc N L)] _1/3 fcNL, 



A£ l (*nl) = /nl[A 2 (/c l )]. 



(5) 
(6) 



The extension to models with Q ^ 1 is straightforward 
in the highly nonlinear regime. The / BL (i) oc x 3 '" 2 scaling 
comes because, once a virialized object is formed, the non- 
linear correlations depend on scale factor a(t) just as the 
background density, £nl oc a 3 , whereas the linear correla- 
tions scale as £ L oc a 2 . If we now allow a density-dependent 
growth-suppression factor into the linear growth law, 



£ L cc [ag(n)] 2 , 

the virialized regime becomes 
U(x)KX 3/2 [g(n)}- 3 . 



(7) 



(8) 



The g(Q) factor may conveniently be taken from the high 
accuracy fitting formula of Carroll, Press & Turner (1992): 



g (Q) = -fi m [niC - n v + (i + On/2)(i + n v /70) 



(9) 



where we have distinguished matter (m) and vacuum (v) 
contributions to the total density parameter. 

The problem is now reduced to one of running a number 
of iV-body simulations to obtain the full form of / NL and 
to investigate its dependence on cosmological model. The 
power of the original HKLM method was the belief that / N l 
was a universal function. Using this assumption, PD fitted a 
^-dependent nonlinear function to a restricted set of iV-body 
simulations - a procedure which works reasonably well for 
spectra with n in the region of —1. However, with the larger 
library of simulations studied here, it is possible to see that 
there is a dependence on the linear spectrum, but one that 
can be described well by a simple fitting formula. 



3 NUMERICAL DATA AND FITS 
3.1 iV-body code 

The carrying out of the required set of numerical experi- 
ments has become much easier recently, partly due to in- 
creased power of computer hardware, but mainly through 
the generous distribution of the Adaptive Particle-Particle 
Particle Mesh (AP3M) code of Couchman (1991). This 
solves the Poisson equation on a mesh to obtain the large- 
scale force, which is then supplemented by exact pairwise 
forces from the near neighbours. Such P3M codes normally 
slow down for highly clustered distributions, but this is 
avoided in Couchman's method by adaptive regridding of 
the densest regions on a finer mesh. 



For the present purposes, the main problem with 
Couchman's code is that it is supplied for an Q — 1 cosmol- 
ogy only. The equations of motion use positions, x, measured 
in units of the box size, Lbox (plus a scaling with the Fourier 
mesh size) . The velocities are defined with these length units 
and a time unit of the age of the universe at the start of the 
simulation. Also, a 'time' variable, p, is used to integrate the 
equations of motion: 



(10) 



In principle, the index a can be altered according to the 
spectrum; we retained a — 3/2, so that both p and the 
scalefactor a were unity at the start of the simulation. 

The standard equation of motion (e.g. Section 14 of 
Peebles 1980) can then with a little effort be cast into the 
form (dashes denoting d/dp) 



V + 



1 2 T' 

- + — + 7T7 

p ap I 



v = Q m (p) 



2ap 



(11) 



where f is the force vector as calculated by the AP3M pro- 
gram, and T = H(p)ti is the product of the physical Hubble 
parameter and the initial time. The required modifications 
to the AP3M equations of motion are therefore to multiply 
the forces on the rhs by ft m {p) and to use the appropriate 
T(p). Note that a non-zero cosmological constant only enters 
through T(p). 

For an Einstein-de Sitter Universe, T = (2/3)a~ 3/2 . In 
general, we have the exact result (Carroll et al. 1992) 



H(a) = Ho^/flvil ~ a- 2 ) + On(a- s - a~ 2 ) + a' 2 . (12) 
There is also the excellent approximation 



H to = l\l-f\-^ S k] /^ : 



(13) 



where / = 0.7fi m — 0.3fi v + 0.3 and Sk is sinh if / < 1, 
otherwise sin. However, this is not needed, since we are only 
interested in T'/T: 

T_ _ J_ » m (2a~ 2 - 3a~ 3 ) + 2a~ 2 n v - 2a~ 2 

T ~ 2ap fi m (a~ 3 — a~ 2 ) + O v (l — a~ 2 ) + a~ 2 ' ( ' 

The only problem with using these formulae to obtain 
T(a) is Couchman's convention that the reference time to 
at which a = 1 is at the start of the simulation, so that 
fi m and Q v in the above formulae are the initial values, not 
those at the desired endpoint of the simulation, designed to 
correspond to the present epoch. To relate the initial and 
final values of f2 m and fl v , we use 



O m (a) = 



fin 



a + Q m (l — a) + Sl v (a 3 — a) ' 



f2 v (a) = 



a 3 ri v 



a + O m (l — a) + fl v (a 3 — a) 



(15) 



(16) 



These formulae work whether a = 1 is regarded as the start 
or the end of the simulation; the important thing is not to 
mix the conventions. 

Lastly, there is the question of initial conditions. The 
AP3M program generates a realization of a displacement 
field corresponding to the given power spectrum. A possible 
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source of confusion, as usual, lies in power-spectrum units. 
Couchman's function pow relates to the notation here via 



pow = 



2tt 2 



[Lbox k] 3 



A 2 (k). 



(17) 



Zeldovich initial conditions are assumed, so that inital dis- 
placement and velocity are proportional. However, in order 
to be properly in the growing mode, the fi = 1 velocities 
should be multiplied by the growth factor O^; 6 (initial) (ef- 
fectively independent of f2 v : Lahav et al. 1991). 



3.2 Simulations and analysis 

There are many degrees of freedom in possible JV-body runs. 
The physical ones are the power spectrum of interest and 
the cosmological model. The size of the simulation box also 
matters, since this must be set so that the fundamental mode 
does not saturate: 



A 2 (27r/L bo x) < 1. 



(18) 



If this condition is violated, the results may not be reliable 
on any scale, owing to the missing power beyond the box 
scale; we used a maximum value of 0.04. It is also necessary 
that A 2 on the initial mesh scale does not exceed unity, so 
that the Zeldovich method used to set up initial conditions 
does not produce excessive shell crossing. For steep spectra, 
this can imply a more restrictive limit on the final box-scale 
amplitude. 

There remain the numerical parameters of the num- 
ber of particles, the total expansion factor, the number of 
timesteps, the force softening and the required force accu- 
racy. We carried out various experiments varying these pa- 
rameters, to see how robust the results were. We generally 
found that even quite simple simulations would reproduce 
the main features of the nonlinear results here, particularly 
the quasilinear regime. We obtained most of our data in a 
standard configuration of N = 80 3 particles, integrated in 
300 timesteps over an expansion factor of 15. A 128 3 Fourier 
mesh was used, with the initial softening set at one cell, 
held constant in proper terms down to a minimum of 0.1 
cell. Some cases were checked in N = 100 3 runs with 600 
timesteps, but these much longer simulations gave essen- 
tially identical results. 

The only factors which influenced the resulting power 
spectra significantly were particle placement and the ques- 
tion of whether to use proper Gaussian realizations in the 
initial conditions. It is common to start iV-body simulations 
by applying the initial displacement field to particles placed 
on a uniform grid, but there are times when this is not desir- 
able. For example, the original grid is often noticeable even 
at late times in void regions. A cosmetic improvement can be 
made by starting from a set of particle positions which are 
irregular but sub-random, such as the 'glass' described by 
Baugh, Gaztanaga & Efstathiou (1995). A simpler alterna- 
tive is to displace each particle randomly within its grid cell, 
which sets up an n = 2 spectrum of initial perturbations in 
addition to those imposed by the initial displacement field. 
With any of these alternatives, some care is needed in ob- 
taining the power spectrum of the density field. 

As in the case of redshift surveys, the particle distribu- 
tion is Fourier analyzed, and shot noise is subtracted from 



the raw power to allow for particle discreteness (see Pea- 
cock & Nicholson 1991). This would be a correct procedure 
at all times if the initial particle positions were also Poisson 
distributed, but the resulting small-scale fluctuations would 
then swamp the desired physical spectrum. The alternative 
starting schemes cure this problem, but introduce discrete- 
ness fluctuations that are initially smaller than Poisson, so 
that subtracting shot noise underestimates the spectrum in 
the early stages of the simulation. Small-scale mixing rapidly 
cures this problem at large k (effectively as soon as the first 
pancakes form), but the overall spectrum will not be correct 
until the large-scale portion has grown enough to exceed the 
shot power 



A 2 hot = [fc/(2^/L box )] 3 . 



(19) 



For N 



particles, the shot power on the box scale 



is A 2 hot = 10~ 4 ' 6 , and this only becomes negligible when 
the box-scale power is ~ 10~ 3 . For n < 0, the shot power 
rises more rapidly with k than the physical power, and so 
in practice larger values of the box-scale power need to be 
used in order that the quasilinear part of the spectrum is not 
affected by shot-noise subtraction. For n = —2, this requires 
box-scale power at the nonlinear limit of 0.04: simulation of 
more negative indices is not feasible without a great increase 
in the number of particles used. 

Even for the analysis of N ~ 10 6 particles, only an 
FFT is practical, in which case the range of k accessible is 
limited by the size of the FFT array. For a 256 3 mesh (our 
normal limit), the Nyquist frequency is only 128 times the 
fundamental, and even here the power spectrum is affected 
by binning and by aliasing (Baugh & Efstathiou 1994). The 
binning correction is, to second order in k, 

A 2 in = A t 2 ruc /(l + fc 2 B 2 /r2), (20) 

where B is the bin size (= simulation box / 256). Aliasing 
effects depend on the slope of the power spectrum; grav- 
itational instability tends to produce an effective index in 
the range —1 to —2, so this is not a problem in practice. 
Normally this size of FFT was sufficient for our purposes, 
but occasionally we wished to extend the results to smaller 
scales. We then divided the data into a set of sub-cubes and 
analyzed each separately, obtaining an estimate of the small- 
scale power spectrum by averaging. Although these samples 
are individually quite nonlinear on their box scale, a few 
comparisons with time-consuming direct Fourier transforms 
of the complete dataset indicated that the overall power 
spectrum could be recovered in this way over a factor of 
300 in k, allowing the maximum range of nonlinearity to 
be probed consistent with the resolution of the simulations. 
It is interesting to note that high accuracy is required in 
the Fourier analysis. A small error in nonlinear power pro- 
duces a small error in the deduced linear power; however, 
because / NL is so steep in the quasilinear regime, this small 
horizontal shift can result in a gross vertical error in / NL . 
More comfortingly, this means that the forward prediction 
of nonlinear power is quite robust with respect to moder- 
ate errors in / NL , and our simple parameterization of the 
numerical results should give very accurate results in most 
circumstances. 

Finally, to create a proper realization of Gaussian ini- 
tial perturbations, each Fourier mode should have a ran- 
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dom phase and a power exponentially distributed about 
the mean. The alternative is to use exactly the expectation 
power, which has the advantage that the large-scale linear 
portion of the power spectrum does not suffer a large scatter 
owing to the limited number of modes. The results on in- 
termediate and small scales seemed identical in either case, 
and so we used 'improper' realizations. 

In this way, we built up a library of results covering 
— 2 < n < and 0.1 < $7 m < 1, considering both open mod- 
els with $7 V = and flat models with £7 V + $7 m = 1. For the 
same cosmologies, we also considered CDM spectra, which 
are parameterized by the shape flh and the normalization 
as (see PD). For a given cosmology, these are degenerate 
degrees of freedom in that a reinterpretation of the simu- 
lation box length will change both fi/i and as- It therefore 
suffices to fix one parameter and vary the other. We con- 
sidered what is roughly the observed shape, Qh = 0.2, in a 
100 Mpc box, with final values of as in the range 0.4 - 
1. For all these simulations, results are available in the form 
of both the final output time and earlier times. To be sure 
that the results were free from artefacts of the initial con- 
ditions, we used only the last factor of 2 expansion, since 
at least a factor 3 expansion is required in order for initial 
transients to die down (Baugh et al. 1995). This gave a set 
of 48 determinations of / NL to be fitted, from 18 distinct 
simulations. 

The next few Figures give a selection of results. Fig. 1 
shows / NL for a variety of power-law spectra with J7 = 1, 
whereas Figs 2 & 3 shows how these change with decreas- 
ing £7, both without and with a cosmological constant. Figs 
2 & 3 illustrate the point made by PD, that / NL steepens 
as we go to lower-density models. Fig. 1 shows that JMW 
were correct in claiming that there was also a significant 
dependence on spectral index, particularly for spectra with 
n < — 1. Such flat spectra have a larger / NL , but / NL also 
appears to steepen for more negative n, whereas JMW sug- 
gested a fit in which the slope of the quasilinear portion was 
independent of n; in retrospect, this steepening can also be 
seen in their data. Interestingly, however, the spectrum de- 
pendence is less extreme for low-density models, as may be 
seen in Fig. 2: for low densities, the quasilinear portion of 
/nl has a similar slope for both n = and n = —1.5. This is 
an important hint for a way of achieving an improved fit to 
the results. What seems to be happening is that the main ef- 
fect of changing the density is to alter the amplitude of the 
virialized /nl(cc) oc x 3 ^ 2 asymptote. JMW suggested that 
this amplitude was also a function of spectrum, being larger 
for n — —2 than for n — by about a factor 2. This could ac- 
count nicely for the different degrees of density-dependent 
steepening: for n ~ —2, the asymptote is at a sufficiently 
high level that raising it still further by lowering the density 
has a small effect on the function. Conversely, for n ~ 0, the 
nonlinear function saturates quite early (for $7 = 1) owing to 
the low level of the asymptote. We now try to find a fitting 
formula to see how well this insight works quantitatively. 

3.3 Fitting formula 

The fitting formula is the same form as used by PD: 




/nl(i) 



1 + B/3x + [Ax]' 



l + ([Ax] a g 3 (Q)/[Vx 1 / 2 ])l 3 



1//3 



(21) 



Figure 1. The generalized HKLM function relating nonlinear 
power to linear power, for an S7 = 1 Einstcin-de Sitter universe 
and power-law spectra with n = 0, — 1, —1.5, —2. The solid lines 
show the fitting formula of Section 3.3. For spectra with n > — 1, 
the function shows very little spectral dependence. For flatter 
spectra, the nonlinear power at a given linear power is higher; in 
particular the slope of the quasilinear portion around A^ L ~ 10 
increases as n decreases, from approximately / m (i) cx x 3 for 
n = to /nl(i) oc x 4 for n = —2. 



This contains 5 free parameters, each of which is potentially 
spectrum dependent. B describes a second-order deviation 
from linear growth; A and a parameterise the power-law 
which dominates the function in the quasilinear regime; V is 
the virialization parameter which gives the amplitude of the 
/nl(x) oc x 3 / 2 asymptote; (3 softens the transition between 
these regimes. 

We proceeded by fitting the Q — 1 results for individual 
power-law spectra, and looking at the trends of the parame- 
ters with n. This suggested a functional form which became 
progressively more nonlinear as n approached —3 (as with 
the power of (1 + n/3) proposed by JMW). The next step 
was to fit what is now a 10-parameter model to the data over 
the range —2 < n < 0, and this was achieved satisfactorily, 
with an overall rms accuracy of about 12 per cent in / NL (j)) 
over the range 0.1 < A„ L < 10 2 5 . Adding data with $7/1 
required very little alteration to the fit, as hoped from the 
discussion in the previous Section: the g term in the fit- 
ting formula seems to be all that is required to incorporate 
different cosmological models. 

The final step was to incorporate CDM results, which 
are important for two reasons. The CDM spectrum is the 
most important example of a spectrum which curves slowly 
so that the effective power-law index n e s = din P/ dink 
varies with scale. It is also a useful case for the present in- 
vestigation, since it contains very little large-scale power, 
but has an effective n that tends to —3 on small scales. 
Pure power-law spectra with n < —2 are very hard to sim- 
ulate: being so fiat, they tend to saturate the fundamental 
mode before enough small-scale nonlinear evolution has oc- 
curred to erase the initial conditions. JMW suggested that 
the nonlinear behaviour of CDM models could be modeled 
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Figure 2. The generalized HKLM function relating nonlinear 
power to linear power, for f2 m = 1 & 0.2 and zero vacuum energy, 
(a) an n = spectrum; (b) n = —1.5. The solid lines show the 
fitting formula of Section 3.3. As U decreases, the nonlinear power 
increases and the quasilinear portion steepens. However, these 
changes are largely associated with the increase of /nl in the 
virialized regime with A^ L > 100. 

via the power-law model corresponding to n e s at the non- 
linear scale. However, it seemed to us that the whole phi- 
losophy of the HKLM method is that the nonlinear power 
at fc N L derives from the linear power at the smaller fc L - One 
would therefore expect that the appropriate treatment for 
CDM models would be to use a different / NL at each k L , ac- 
cording to the tangent spectral index at that point. This is 
an important assumption, because it means that the small- 
scale power in CDM models should be representative of the 
n < —2 spectra which are otherwise so hard to treat. With- 
out this assumption, we found our CDM results hard to fit: 
they reach nonlinear powers on the smallest scales greater 
than would be expected from even the n = — 2 power-law 



o 




Figure 3. As for Fig. 2, but for spatially fiat models with 
f!v + f2 m = 1. The generalized HKLM function relating nonlinear 
power to linear power, for Q m = 1 & 0.2. (a) an n = spectrum; 
(b) n = —1.5. The solid lines show the fitting formula of Section 
3.3. The Q dependence is weaker than for open models, consis- 
tent with the idea that all that matters is the growth-suppression 
factor g(£l). 

fits (see Fig. 4), but the nonlinear response at larger scales 
is less extreme, as would be expected if the effective n was 
larger. However, although the trend of nonlinear response 
with scale is as expected, the CDM results generally lay be- 
low those predicted from the power-law fits using n e ff(fc L ). 
This is not unreasonable: if we compare a CDM spectrum 
with its tangent power-law spectrum, the power-law spec- 
trum has a greater amount of power integrated from k = 
to k—k L . A practical means for accounting for this difference 
is to conjecture that the nonlinear behaviour of the CDM 
spectrum will be characteristic of its tangent index at some 
slightly smaller scale, and a shift of a factor 2 in k gives 
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Figure 4. As for Fig. 2, but for standard CDM spectra. The 
generalized HKLM function relating nonlinear power to linear 
power, for (a) for an Q = 1 Einstein-dc Sitter universe; Qh = 0.2 
and ag = 0.3 & 0.5. (b) for an Q m = 0.1 universe with Q„ = & 
0.9; Qh = 0.2, erg = 0.7. The solid lines show the fitting formula of 
Section 3.3. Note that the CDM curves arc steeper even than the 
n = —2 power-law results, suggesting that the small-scale CDM 
behaviour is characteristic of the tangent spectral index there, 
which can be as negative as n = —2.5 for the models studied here. 
Figure 4b shows clearly the extra small-scale power produced in 
the case of an open model by comparison with a a flat model of 
the same Q m . 

outstandingly good results; 

n L (fc L ) = ^(fc = fc L /2). (22) 

The exact size of the shift is not critical, and we have made 
no attempt to treat it as an additional parameter to be opti- 
mized. Very good predictions of nonlinear CDM spectra are 
achieved even simply using the unshifted tangent spectral 
index. A shift in this sense is likely to be required for any 



convex spectrum, where n e s decreases as k increases. Since 
this prescription clearly also works for pure power-law spec- 
tra, we propose this as a general method for dealing with 
any smoothly-curving convex spectra that are hierarchical 
in the sense that A 2 increases as k increases. The obvious 
exceptions are therefore spectra with a small-scale cutoff, 
such as hot or warm dark matter, and these will require 
separate treatment. 

We therefore performed a global fit to the power-law 
plus CDM data with this assumption, and obtained an ex- 
cellent fit over all cosmological models, with an rms accu- 
racy of about 14 per cent in / NL (x). This a demanding level 
of agreement, since / BL (i) is so steep; the scatter in the 
transverse direction, which governs the accuracy of power- 
spectrum reconstruction, is only about 7 per cent. As dis- 
cussed above, errors in / NL partially normalized themselves 
away when predicting nonlinear power, and this figure of 7 
per cent is also the approximate rms accuracy with which 
this method here will predict A^ L for a given linear spec- 
trum. The best-fit parameters are 



0.482 (l + n/3)" ' 947 


(23) 


0.226 (l + n/3)" 1 - 778 


(24) 


3.310 (l + n/3)" - 244 


(25) 


0.862 (l + n/3) -0 ' 287 


(26) 


11.55 (l + n/3)- ' 423 . 


(27) 


Note once again that the cosmoloj 


peal model does not 



enter anywhere in these parameters. It is present in the fit- 
ting formula only through the growth factor g, which governs 
the amplitude of the virialized portion of the spectrum. This 
says that all the quasilinear features of the power spectrum 
are independent of the cosmological model, and only know 
about the overall level of power. This is not surprising to the 
extent that quasilinear evolution is well described by the Zel- 
dovich approximation, in which the final positions of parti- 
cles are obtained by extrapolating their initial displacements 
by some universal time-dependent factor. All information on 
the cosmological model is hidden in this extrapolation fac- 
tor, and therefore the model should have no effect if we scale 
to displacements of the same size. The power spectrum in the 
Zeldovich approximation has been calculated analytically by 
Taylor (1993) and by Schneider & Bartelmann (1995), and 
it would be of interest to compare their results with ours. 



4 DISCUSSION 

We have investigated in detail the scaling formula of HKLM 
for the evolution of clustering statistics in cosmology. Al- 
though not completely spectrum-independent, their ap- 
proach can be made to give a good fit for a variety of spec- 
tra, provided one uses a simple dependence on the tangent 
slope of the linear power spectrum. Because of the need for 
a spectrum-dependent correction, we have not provided a 
fitting formula for £, nor for the inverse nonlinear function. 
Inverting observed nonlinear data is now in any case an it- 
erative process, and £ is a statistic of less practical interest 
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Figure 5. The present fitting formula for the nonlinear function 
(solid lines) compared with the spectrum-independent form sug- 
gested by PD (dashed lines), for the cases n = —1 & —1.5 and 
open models with Q m = 1,0.5,0.2. The PD formula was an ap- 
proximation to the average effect of these different spectra, but 
the detailed agreement is often poor, particularly at low f2 for 
flatter spectra. However, these deviations are small in the regime 
of the data used by PD (A& L < 3). 

than the power spectrum. It is probably best to proceed nu- 
merically with the forward nonlinear function as a starting 
point. 

The main feature of the nonlinear HKLM function in 
this work is a power-law on intermediate scales which is 
rather steep, /nl(:e) oc x 1+a with a ~ 3.5 - 4.5, and it is 
a challenge to understand this result. Padmanabhan (1996) 
has given arguments to suggest that the intermediate slope 
should be /nl(x) oc x , but it seems that this is not the 
true value. For spectra with n ~ 0, such an index works well 
enough, but this appears to be an artefact of the relatively 
rapid onset of the virialized regime. For lower densities or 




very negative n, the virialized regime occurs at larger pow- 
ers, so that the steeper intermediate behaviour of / NL (a;) 
is seen more clearly. This steep quasilinear function has an 
interesting implication for the scaling of small-scale power 
spectra and their evolution with time. For A 2 2> 1, we have 



&NL — [A NL ] ^ 

^nl oc [D 2 (a) Al] 1+q \ 



(28) 



where D(a) is the linear growth law for density perturba- 
tions. For a power-law linear spectrum, this predicts a quasi- 
linear power law 



AL ocZ?( 6 - 2 «( 1+ «)/ 3 e 



(29) 



where the nonlinear power-law index depends as follows on 
the slope of the linear spectrum: 



0-- 



3(3 + n)(l + a) 
3 + (3 + n)(l + a)' 



(30) 



For the observed index of f3 — 1.8, this would require n ~ 
—2.2, very different from the n = that would give f) = 1.8 
in the virialized regime. However, especially for low density 
models, the virialized regime is only reached on very small 
scales and the observed clustering data are dominated by 
quasilinear effects. It is then interesting to note that the 
predicted evolution is faster than the linear A oc D (a); 
this may be of relevance in understanding the weak angular 
clustering of faint galaxies (e.g. Efstathiou et al. 1991; Roche 
et al. 1993). 

These results should be of practical use for a variety 
of cosmological investigations. The most obvious case is the 
theme pursued in PD: linearizing observed clustering data 
in an attempt to infer the underlying linear power spec- 
trum. The fitting formula used by PD was restricted by the 
assumption that / NL was spectrum independent, and the 
results given here are of considerably greater accuracy, par- 
ticularly for n < — 1. How much this matters depends on the 
application of interest; because / NL is so steep, it is possible 
to predict A 2 L rather badly and yet be able to infer A 2 with 
tolerable accuracy. This is illustrated in Fig. 5, which com- 
pares the present fitting formula with the one given in PD. 
There are deviations of up to a power of ten in A^ L in the 
case of great nonlinearity and low density. However, in the 
regime of the data actually used by PD (A 2 ^ < 3), the errors 
in the deduced A 2 from a given nonlinear power are at most 
around 20 per cent; the work reported here thus does not 
imply any serious revision of the conclusions reached in PD. 
Nevertheless, there remains the challenge of understanding 
the highly nonlinear portion of the the power spectrum, and 
the improved fitting formula from this paper should be of use 
in attempts to interpret the clustering data in this regime. 



ACKNOWLEDGEMENTS 

SJD was supported by a SERC/PPARC research stu- 
dentship during part of this work. We salute Hugh Couch- 
man for his tremendous contribution to cosmology in making 
his AP3M code freely available. We thank Bhuvnesh Jain for 
helpful correspondence on this subject. 



8 J. A. Peacock and S.J. Dodds 



REFERENCES 

Baugh CM., Efstathiou G., 1994, MNRAS, 270, 183 
Baugh CM., Gaztanaga E., Efstathiou C, 1995, MNRAS, 274, 
1049 

Carroll S.M., Press W.H., Turner E.L., 1992, ARA&A, 30, 499 
Couchman H.M.P., 1991, ApJ, 368, L23 

Efstathiou C, Bernstein C, Katz N., Tyson T., Guhathakurta 

P., 1991, ApJ, 380, 47 
Hamilton A.J.S., Kumar P., Lu E., Matthews A., 1991, ApJ, 374, 

LI (HKLM) 

Jain B., Mo H.J., White S.D.M., 1995, MNRAS 276, L25 (JMW) 
Lahav O., Lilje P.B., Primack J.R., Rees, M.J., 1991, MNRAS, 
251, 128 

Padmanabhan T., 1996, MNRAS, 278, L29 
Peacock J.A., Nicholson D., 1991, MNRAS, 253, 307 
Peacock J. A., Dodds S.J., 1994, MNRAS, 267, 1020 (PD) 
Peebles P.J.E., 1980, The Large-Scale Structure of the Universe. 

Princeton Univ. Press, Princeton, NJ 
Roche N., Shanks T., Metcalfe N., Fong R., 1993, MNRAS, 263, 

360 

Schneider P., Bartelmann M., 1995, MNRAS, 273, 475 

Taylor A.N., 1993, in Proc. Cosmic Velocity Fields, 9 th IAU Conf, 

cds F. Bouchet, M. Lachieze-Rcy, Editions Frontieres, Gif- 

sur-Yvette, p585 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science TjgX macros. 



